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1. Introduction 

The AdS/CFT correspondence Q conjectures an equivalence between string theory back- 
grounds and gauge field theories. The original duality proposes a correspondence between 
the 3 + 1 dimensional N = A supersymmetric Yang - Mills (SYM) theory and type IIB 
string theory on AdS^ x . The conjecture is extended by further dualities between differ- 
ent string backgrounds and field theories with different properties than A^ = 4 SYM, such 
as spacetime dimensionality, less number of supersymmetries, etc. The string background 
studied in this work is a generalisation of the one proposed in |^ and utilised in Q for the 
construction of a string dual to N = 1 supersymmetric QCD-like theories. 

The string duals of minimally supersymmetric QCD-like field theories employs D5- 
branes with world-volumes extended along the Minkowski directions and wrapped on a 
compact two - sphere. Matter is added by the fiavour D5-branes separated from the 
colour branes along a radial direction Then the theory includes matter fields in the 
fundamental representation, which is more suitable to phenomenological studies. Even 
though the exact solutions to this background are so far unknown, solutions for specific 
values of Nf and Nc-, possible asymptotics and numerical solutions are already available and 
their properties have been studied to some extent. For some papers where the fiavoured 
MN background has been studied see ^, 0, ^, ^ • For other papers dealing with flavours 
added by backreacting fiavour branes see [0, |ll|, |l|, ||, |l7|, . 
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One of the known solutions has the property that the couphng constant in its dual field 
theory has the so called walking feature |21]. This is a required feature for the coupling 
between the techniquarks of technicolor models that are used as a natural model for the 
electroweak symmetry breaking. Similar to the case of QCD, the coupling constant diverges 
at the low energies. In this work a new numerical solution to the background with a walking 
gauge coupling is presented. The major difference to the one in |^ is that the coupling 
constant remains bounded everywhere. This is achieved by adjusting some parameters so 
that the background functions have a certain asymptotical behaviour. These asymptotics 
have already been classified in |Q H], [2^ . 

The second section is a review of the background with its known solutions. The prop- 
erties of the known solutions is presented as well as the classification of the asymptotics. In 
the third section, the bounded walking solutions are constructed by imposing the solutions 
to have certain classes of asymptotics. It turns out to be that there is a constraint on the 
parameters of the system to produce the desired kind of solutions. In the fourth section, 
the Wilson loops have been calculated for the presented solution. The heavy quark po- 
tential exhibits a phase transition. Similar phase transitions of the heavy quark potential 



have been related to the phase transitions of the van der Waals gases [^] and [16[. Finally 
in the fifth section, not relevant to the bounded walking solutions, the effect of flavours is 
studied by finding flavoured solutions in a perturbative way. 

2. Review of the string dual to J\f = 1 SQCD-like theories 
2.1 Description of the background 

In this chapter the background that is proposed to be dual to SQCD like theories will 
be reviewed very briefly. It is built on the Maldacena - Nuiiez background, which is a 
type IIB supergravity solution introduced and proposed to be dual to a = 1 SYM-like 
theory in 0] ( see also |24| where a 4d solution was found). This theory includes Nc colour 



D5-branes extending along the Minkowski directions xi^s and the compact directions {6, (p) 
that provide the colour symmetry of the theory. Then, Nf flavour D5-branes smeared along 
four compact coordinates [9, if, 6, ip) are added to introduce an open string sector. Their 
worldvolume extends along the Minkowski directions, the radial coordinate that will be 
denoted as p and a compact direction ^/J. The open string sector gives rise to fundamental 
fields in the field theory. This system is governed by the action S = Sub + S flavour which 
consists of the Type IIB supergravity action and the Dirac - Born - Infeld + Wess - 
Zumino (DBI-I-WZ) action for the flavour branes. They read: 



SiiB =:r^ I d^'x^ 

^'^(10) 

Nf 

''flavour 



V J Me ^ J Ma J 



(2.1) 



where is the pullback of the metric and ^[Ce] s the pullback on the worldvolume of 
the RR six-form of the background. Having added the flavour branes, their effect on the 
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background cannot be neglected, thus they are said to be backreacting. The action of the 
D5 flavour branes is six dimensional and they live on constant values of the two spheres' 
coordinates {9,(j)) and {0,4)). The smearing mentioned in the above paragraph avoids the 
need of delta functions for these six dimensional objects in the ten dimensional theory and 



this technique was first used in |18|. The ansatz for the metric of the string background 
is a generalisation of the background in [Q]. The metric describes a space of the topology 
M}'^ X M X 5^ X 5^ and it reads in the notation of ||2H: 



dx\^ + Y{p) {Up'^ + (W3 + Co^f) + ip(/>) sinh (t(/>)) (witJi - UJ2^2) 
+ i (P(p) cosh (t(p)) + Q{p)) {loI + u^l) + \ (Pip) cosh (r(p)) - Q(p)) {CjI + 

(2.2 



The dilaton 4>{p) is a function of the radial coordinate and the RR 3-form is given by 

Nr 

-^(3) =^ [- i^i + Kp)dG) A (cD2 - Hp) sin {e)dip) A (ws + cos {e)dip) 

+b'{p)dp A {-de A lDi + sin {9)d(f A UJ2) + (l - b'^{p)) sin {9)d9 A dip A 0)3] (2.3) 

Nf / ^ 

- sin {Q)d9 A dip A {dip + cos {9)dLp 

where 

a)i = cos {Tp)d9 + sin (-0) sin {9)d(p, 

UJ2 = — sin {'ip)d9 + cos (V') sin {9)d(p, (2-4) 
u)i = (i^/; + cos {9)d(p. 

The background functions P, Q, 1", 0, a and h are determined by a set of equations. These 
are BPS equations which are first order differential equations obtained from the super- 
symmetry transformations of the fermions and they ensure that the background satisfies 
the equations of motion of the gravity plus flavour branes action. After making some 
redefinitions like a = -5— tt-^-Vt: b = -?r and the introduction of the function 

( 2Nc-N, \ 

u; := (J - tanh(r) Q + — (2.5) 
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the BPS equations are written as follows |22]: 



P' = 8Y -N 



d. 



Q 



^ ^ cosh(r^ 

On 



r'(/9) + 2sinh(T) 



(2iVc -Nf) 2uj , ^ o, 
2 cosh(r) 

ley 

P2 _Q2 

2Q cosli(r^ 
2a; 

P2 cosh(T) 



(2.6) 



P"sinh^(r) + Q Q + 



2N, - Nf 



with an algebraic constraint u; = and ^> = (P^ — Q^) Y'^e^'^. One then calculates r and 
Q using (|2.6|) as: 



cosh(r) = coth(2p) 



Q=[Qo + 



2Nr - N 



cosh(r) + 



\ ^ (2pcosh(r)-l) 



(2.7) 



These equations can be cast into a single differential equation for P that reads |£2[: 

'P' + Q' + 2iV/ P'-Q' + Nf 



P" + (P' + Nf) 



+ 



4cosh(r) = 



(2. 



P-Q P + Q 

Once having obtained P, the other functions can be calculated in terms of it using the 
equations (2^). The equation ( |2.8| ) will be referred as the master equation. Its detailed 
derivation can be found in [BTl 1231 . 



2.2 Type A backgrounds 

In particular, the backgrounds described above represent a relatively general case and 
they are called type N backgrounds in the notation of Q. There is a special case in 
which the backgrounds are called type A. In this case, the fibration of the S'^ and the 
becomes simpler. On the field theory side, as opposed to type N backgrounds, the gaugino 
condensate in the field theory vanishes. It is possible to reproduce the type A backgrounds 
by setting a=T=0. This unifying approach was presented in |22|. 



When dealing with type A backgrounds, it is useful to redefine the background func- 
tions: 
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In this notation the metric and the 3-form of the type A backgrounds are given by: 
dxl 3 + 4Y{p)dp'^ + H{p){de^ + sin^{e)dip'^) + G{p){de^ + sin^Od^'^^ 



ds^ = e<^(^) 



+Y{p){dij + cos{9)d(f) + cos{e)dipf 



(2.10) 



^(3) 



-d [a{p){uJi A (Di — W2 A 0)2] 



Nf-Nc Nc , , X . X 
uji AL02 + — cui A (I'2 ) A (0^3 A 0)3) (2.11) 



with LUi = d6, L02 = sm{9)dip, = cos{6)dip. The master equation in type A backgrounds 
reduces to 



H" 



ldpH + l{Nf-N, 



d,H + N 



f 



Nf + 2dpH 



H 



H H — — p - C 



0. (2.12) 



As before, the other type A background functions can be calculated in terms of H{p) and 
integration constants 0]. 



G = H + 



Nf - 2N, 



p-C 



2 

y = \d,H + -(Nf-N,) 



^-0 + 



Nf-N, ^ iV, 



AH 



AG 



dp 



(2.13a) 
(2.13b) 
(2.13c) 



2.3 Some known solutions 



Although the master equation is a highly non-linear second order differential equation, it 
permits simplifications by setting its parameters to certain values. For these specific cases 
there are some known analytic solutions. The ones that are considered in this work are 
listed below. 

While only P{p) or H{p) is given here, it is possible to calculate the remaining functions 
by the equations and the definitions of the previous section. For type N backgrounds see 
the equations (^]^ - 2X). Equations (|2.13|) are valid for the type A case. It is always 
possible to translate between the "P-Q" and notations using equations (|2.9D . 



Unflavoured solutions Nf = 
• A type N solution from is: 

P = 2Nc{p - po), Qo = -N, - 2N,po, p>po>-(X) 



(2.14) 



The type A limit (r, cr ^ 0) of the type N solution in eq. (I2II ) above IS 



H = Nc{p - p*), p* : = 



2A^. 



Nr 



, P> P* 



(2.15) 
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Solutions for Nf = 2Nc 

Nc = 2Nf is a special case of the background where the functions have some different 
properties. In this case Q{p) does not diverge unhke the other cases but it goes to a 
constant. 

• A type A solution called the "conformal" solution in IH reads 

(2.16) 

where ^ is a real number between and 4. 

• A deformation of the above solution 

^ = ^(9±3) + ^e^ (2.17) 

also solves the master equation and it is of type A. 
Further solutions for this case have been presented in 

An arbitrarily flavoured type A solution 

The following type A solution is an obvious solution of equation ( 2.12| ). However it is not 



physical since the background function Y(p) turns out to be zero (seeeq. (^l3D). 

Nf - Nr 

H{p) = ^^{c,-p) (2.18) 
2.4 IR and UV asymptotics and their classification 

Besides the exact solutions that are valid for specific cases, there are asymptotic expansions 
written for the generic cases. There are two kinds of asympotics for the UV and three for 
the IR limit. The UV asymptotics are referred as "Class I and 11" while the IR ones are 
"Type I, II and III". 

UV asymptotics 

As discussed in |22], equation (3.17) of |22| implies that the P must diverge as p — > oo 



4£ 

unless Nf = 2Nc. There are two possibilities. P either asymptotes to oc p or oc e a . They 
are called Class I and II UV behaviours, respectively. 

• Class I In cases where Nj ^ Nc, the class I asymptotics are of two kinds depending 
on whether Nf is larger or smaller than 2A''c. For Nf > 2A''c, H{p) asymptotes to 
\iNf — Nc), while for the other case it asymptotes to ^(2A''c — Nf)p. This class of 
asymptotics contains terms up to the first order in p. The leading term of G and Y 
as p ^ oo is a constant, namely The dilation like e^*^ oc e^P / p. 

• Class II In class II solutions, H{p) and G{p) asymptote to c+e"3" regardless of the 
parameters. G, Y are also proportional to e~ and the exponential of the dilaton e^"^ 
approaches a constant. 
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Since r vanishes at the large values of p , there is no distinction between type N and type 
A backgrounds in the UV region. 



IR asymptotics 

The three types of IR behaviours are classified by the vanishing of some of the background 
functions in the IR region. The definition range of p, thus its smallest value called the IR 
value varies from case to case. The asymptotics of special type A backgrounds are reviewed 
here, since it is the relevant case for this work. In type II backgrounds the IR limit is either 
p = or p = — oo. 



Type I: In such backgrounds the IR limit is defined as /> — > — oo. The series expansion 
for H has the following form: 

H{p) = ^^^y^(ci - p) +Y.V,{p)e'^^ (2.19) 

k>l 

where Vt are polynomials of order A; + 1 in p. From this Y is found to be proportional 
to e 3 . As stated above, the solution given in eq. (|2.18D is a bad solution because of 
the vanishing of Y{p) Vp . The type I solutions are solutions that do not suffer from 
this, however they approach eq. (|2.18| ) as p ^ oo. 

Type II: In backgrounds with type II IR behaviour, the IR limit is defined as p — > 0. 
The series expansion for the background functions H, G, Y, cj) is either of the forms 



^ ^ A^f-A^c 1 3C{Nf - Ncf - h'^{2Nc + Nf) ^ ^ ^ 

(p = <Po+ V ^+^J UnJ ^-P+- forC<0 

2hi llChf 



or 



^^C + ..p.-^| + fjp + ... 

G = /.,P^-(^ + ^jp+--- 

hi I (hi NA , 72C2 + 20/if + 3C(10Ac - 3A/) i 



4p 



Nc 1 3C - hl{3Nf - 2Nc) 
2hl^' ^ 12Chl 



^o + ^P^+ ^' P + --- forC<0 



(2.20) 
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depending on the sign of C. Note that one of the functions vanishes at /) = while 
the other is non-zero, C. For C = 0, both of the functions are zero and these kind 
of backgrounds are called type III. Their leading behaviour is like pa as described 
below. 

• Type III: Finally, there is the third kind of IR asymptotics that read. 

H = hip-^ +^{5N,-7Nf)p+•■■ 
^^ (2.21) 
G = hip^ - —{5N, + 2Nf)p + --- 

Type III asymptotics are valid instead of type II when C = 0. Both H and G are zero in 
the IR limit. 

The above review summarises the known possible asymptotical IR and UV behaviours 
of the solutions for the type A BPS equations of the background being studied. A new 
numerical solution found with class II UV and type II IR behaviour as well as the conditions 
for its existence and some of its implications on the field theory side such as the gauge 
coupling will be presented in the following section. 



3. Bounded walking solutions of the background 

Walking is a desired property of the so-called techniquark coupling constants of technicolor 
models that are used as a natural model for the electroweak symmetry breaking. Because 
of phenomenological reasons their running coupling constant must differ from the QCD 
coupling constant. The techniquarks need to remain strongly coupled for some higher 
energy scale. In other words, the coupling constant is needed to "walk" at a certain value 
before it decays to zero in the far UV. This kind of walking has already been observed in 



the background studied in this paper |21]. The coupling constant can be translated for 



type A backgrounds from the background functions as (see 

1 



While the coupling constant in |21| diverges in the IR like the QCD coupling, it can be 
useful to find backgrounds with bounded field theory coupling constants, in backgrounds 
with type II IR and class II UV asymptotics is expected to be constant for a certain interval 
and drop quickly at a certain point in the UV. That is indeed what happens at least for 
some values of the parameters hi and C in the IR expansion. 

To investigate the parametrical constraints to obtain the desired kind of solutions a 
set of sample parameters 

(C,/ii) G {1,10,50,500,5000} x {0.01,0.7,0.78545,0.78552,0.78657,0.85,10,15} 

has been considered where Nf has been taken as zero at the first step. These solution is 
obtainable with some of the combinations of parameters from this set. Figure |l| contains 
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H[p] <3[p] *[P] 




Figure 1: A sample walking solution of the background with type II IR and Class II UV asymp- 
totics. The plots are for C = 5000, hi — 0.7852. There is no significant change in the coupling 
constant until the far UV values of p. 

the plots of the background functions and the coupling constant for such a solution. The 
functions H(p) and G{p) have a square root behaviour in the IR while they diverge expo- 
nentially in the UV region. The dilation has no singularities and asymptotes to a constant 
in the UV. The coupling constant has the desired shape. The table in figure § contains plots 
of the coupling constant for different parameters C and hi. As it can be seen from that 
table, the existence of the desired solutions depends on both of the integration constants. 
Moreover, the bound of one parameter is related to the value of the other. 

3.1 Parameter constraints for the desired kind of solutions 

A more precise list of observations that can be concluded from figure ^ is the following: 

• For solutions of common hi, whose coupling constants blow at a certain point in the 
UV, become a smoothly and rapidly decaying behaviour as C increases. At even 

4p 

larger C, they turn into Class II solutions (oc e 3 ). See the change of the yellow 
{hi = 0.78552) solution for different C. 

• For fixed C, smaller hi is required to have longer plateau. However, there is a 
minimum for hi that decreases as C increases. 

• While for the {hi = 0.85) and {hi = 0.78657) lines C = 5000 results in a longer 
plateau than C = 500, for the solutions {hi = 0.78552) the case is opposite. 

Adding Nj 

To get an intuition on the effect of flavours, a similar investigation has been made for 
flavoured solutions. Figure |3| contains plots of solutions with two different C and constant 
hi = 0.78659. Different lines correspond to different x = ^ £ {0,0.5,1}. Cases with 
Nf > are only nice for larger C than the cases with Nf = 0. However, once C is large 
enough the value of x has no effect on the solutions. 
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Figure 2: The plots of the field theory coupling constant for diff'erent values of C and hi Each 
cell contains plots with common C . The colours represent the value of hi 



4. Wilson loops 

Wilson loops are gauge invariant operators that allow to compute the energy of an quark 
- anti-quark pair. Their relation to gauge / gravity conjecture has been stated in |^ and 
the calculation is reviewed by |2G] in detail. The Wilson Loop operator constructed over 



a closed space-time loop C representing the world-line of a quark - anti-quark pair. The 
world-line of pair we consider forms a rectangle where the pair is created in one side of it 
and annihilated in the other. When T, the time leg of the rectangle we consider is very 
large, the energy of the pair as a function of the separation of the pair can be estimated 
as follows: 

< W{C) >= A(L)e-^^(^) (4.1) 



On the other hand, it is proposed in |^5[ that the expectation value of the Wilson loop 
operator can be calculated from the Nambu - Goto action for a string stretching along the 
radial coordinate p. The endpoints of the string trace the quarks' world-lines living at the 
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Figure 3: On the left hand side are the plots for C = 500, the flavoured solutions do not look as 
desired. The right hand side plot represents the solutions with C = 5000 however, one sees that all 
of them are good and exactly the same. 



UV region (large p). 

< W{C) >oc e"^ (4.2) 

I dTda^\det{Kp)\. (4.3) 

where hap = g^udaX^dijX^ is the induced metric on the string world-sheet. After fixing 
the world sheet parametrisation tor = t and a = x), the determinant of the induced metric 
for the static string stretched to along the p - x plane is: 

det{hap) = gttgxx + gttQpp {Oxpf (4.4) 
= a'^e^-^ (-l-4yp'2) 

where p' = Since there is no t dependence in the integrant, the t integration can be 
performed straight away to obtain: 

S=^ [ dxe'^yrTiy^ = I_ [ (4.5) 
2tt J 2iT J 

The Lagrangian does not depend explicitly on x, therefore the x - Hamiltonian 

n = f/-C (4.6) 

is a constant of motion. In fact, it is equal to — e'^^''"-' where po is the point on which p' 
and consequently ^ vanishes. Using this it is possible to write 

dp 1 Ve2<^W + e2<^(^o) 

— = — (4 7) 

dx 2i/y(p) e't'^Po) ^ ' 

to express the quark separation in terms of pQ, the deepest coordinate that the string 
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reaches: 



pi , g</'(po) 



L= dx = 2 2y^)^=====dp. (4.8) 

Similarly, writing the energy of the pair from the equations ([4.1|), ( |4.2| ) as a function of pQ 
as 

/■pi J™ /"Pi 
y= / C^dp-2x / 2yy(^e<^('')dp (4.9) 
JpQ dp Jq 

one obtains a parametric relation between the quark - anti-quark separation and energy. 
In equation ( [4.9D , two quark masses calculated as the length of a string stretching from the 
bottom of the space to the quarks, have been subtracted from the energy to consider just 
the potential energy of the pair. 

Here we calculate the Wilson loops on a background with bounded walking coupling 
constants described previously. We take the loop to be the on the x - t plane at a very 
large radial coordinate p = pi. Figure § shows the plots of energy - separation relation 
for configurations with quark masses mg « 500 and tuq w 3700. One sees two branches 
and the lower one has linear confinement. The upper non-linear branch is less stable due 
to its larger energy. The branches coincide at a certain point at a very high po ^ Pi after 
which the string cannot stretch any further. This is interpreted as string breaking. Even 
though there are no flavour branes in this background to allow formation of dynamical 
quarks where the string breaks, the derivative of the background function H{p) is singular 
at the origin. That affects the energy and the separation integrals and thus is said to allow 
phenomena as if there were flavour branes. 

Such behaviours have been observed and discussed in similar solutions of this back- 
ground, see for example 23, 27]. In [^], [27| and [^] flavoured systems with massive 
dynamical quarks are found out to have the energy - separation relation that qualitatively 
resembles the phase transition in the Gibbs free energy - pressure relation of van der Waals 
gases that involves "quasi-static phases" . Therefore, it has been argued that it can be 
interesting to study quark potentials in gauge theories with string duals by relating them 
to the phase transitions of van der Waals gases. 

In the present background, another interesting phase transition occurs if one puts the 
brane where the heavy quarks live on at distances of smaller values of p. In other words, 
one reduces the mass of the quarks. Then a phase transition occurs at the configurations 
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Figure 5: Decreasing tlie lieavy quark masses, tire phase transition at p < pi is replaced by another 
transition at a lower value of p. The plot to the left shows the general picture while a zoom of the 
transition region is presented on the right hand side. 



MQ=224 .139 



-376 . 02 
-375.03 
-376.04 
-375.05 
-376.06 
-376.07 




-800 

-820 

-840 „ 

-860 

-830 

-900 



-\ PO 



-2000 
-4000 
-6000 _ 



1 2 3 4 5 



Figure 6: The relation between the quark potential energy V and the deepest coordinate that the 
string reaches po- The plot to the left belongs to the case with Mq k, 3700, the plot in the centre 
belongs to the case with Mq « 500 while the one to the right belongs to the case with Mq « 225 



with about p O.lpi. Figure ^ contains the plot of the entire picture and the zoom of the 
region of transition 

Figure ^ shows the relation between the quark potential and the deepest coordinate 
that the string reaches for the three cases mentioned previously in this chapter. The phase 
transitions occur at the points where the potential reaches a maximum or a minimum. 

Concluding this chapter, it is possible to say that one observes phase transitions of the 
heavy quark potential energy also in the bounded walking solutions of the background we 
presented in this work. Where the transition occur and what their qualitative properties 
are, depends on where one places the heavy quarks. None of the phase transitions here is 
exactly of the same type as the one of the van der Waals gases. However, there is still a 
possible relation since the pattern matches partially. 

5. Perturbative addition of flavours 

For the non-flavoured cases there are exact solutions known to the background functions. 
The idea of this section is to look for small x = expansions of some flavoured solutions, 
whose X ^ limit are the known non-flavoured solutions. In other words, we perturb the 
non-flavoured solutions around x = and plug the expression 

H{p) = Ho{p) + hi{p)x + h2{p)x^ + ... (5.1) 

into the H - equation and solve the differential equations obtained for hn{p) functions. The 
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observation is that the properties of the solutions change with the addition of flavours. 



Backgrounds with H{p) = ^ 



As mentioned around equation ( 2.1^ ), this solutions are known to be unphysical Q since 



they have vanishing background function Y = Q for all p and therefore cannot be considered 
as a 10 dimensional background. However, as soon as one adds flavour corrections up to 
the /cth order in x, one obtains correction terms of the form 

where Vk are polynomials of order A; + 1. This solutions are as reviewed before known 
as type I and they are well-behaved solutions. Thus we say, a solution that is not even 
a background turns into a physical background after the addition of flavours. That the 
solution H = p/2 being the limit of type I solutions as /? — > — oo has already been discussed 
in i. 

Backgrounds vi^ith H{p) = p 

Another solution to the H-equation is H{p) = p. From this solution one calculates a dilaton 
that diverges at p = 0. Diverging dilatons in the IR are considered as bad singularities 



causing unphysical effects - see section 5 of . That is because e"^ has to be finite in the 



IR limit. Following the perturbative flavouring procedure described above, one considers 
correction terms that come in powers of x and solves the differential equations for the 
coefficients that are functions of p. Then one obtains the solutions of the form known as 
type II. In type II solutions the dilaton does not have a bad singularity [^]. An example 
for a type II dilaton is depicted in figure ||. 

Backgrounds vi^ith P = 2p 

This unflavoured type N solution is known to be a smooth solution that has no bad sin- 
gularities. Perturbing this solution with flavours brings it to the following form (up to the 
first order in x and second order in p). 



P{p) = 2p + xyj^ + cip-- + 0{p') j + 0{x') (5.2) 
As usual, the dilaton is calculated from the equations ( |2.6D as: 

""'^ (p2_Q2)ysinh2(r) ^^-^^ 

One sees that the dilaton diverges to — oo in the IR after adding the flavours. This diver- 
gence is not considered as bad since e"^ remains bounded. Figure contains the plots of 
the background functions P and the dilaton before and after the perturbative addition 
of flavours. 



-14- 



0(P) 



P(p)=2p, P (p) =2p+pl (p)x 

. 008r 



0.006 
.004 
. 002 




-10 



. 0002 . 0004 . 0006 0^&(TO8 O.OOl 



Figure 7: The change of the background P{p) — 2p with flavours: The plot to the left shows P 
and the one to the right shows 4>. The solid line is the unflavoured solution while the dashed line 
is the solution perturbed with solutions 



6. Conclusions 

In this work a solution to the string background has been presented which is dual to SQCD- 
like field theory with a gauge coupling that has a walking property. It has been observed 
that the integration constants that come from the BPS equations must be be constrained 
by a relation between each other in order to obtain such kind of backgrounds with walking 
gauge couplings. After obtaining the exact numerical solutions it is possible to make 
calculations for the quantities of the field theory such as the quark - anti-quark potential 
using the Wilson loops. The Wilson loop calculations pointed out that the heavy quark 
potential for the field theory experiences phase transitions. The qualitative properties of 
those phase transitions depend on the quark masses which depends on the radial coordinate 
of the D-brane on which the quarks live. Possibly these phase transitions can be linked 
to the phase transitions of van der Waals gases by analogy. Finally, in addition to the 
calculations related to the presented solution, the effect of the presence flavours has been 
investigated. Some already known unflavoured exact solutions have been flavoured by 

Nf 

sl small X = perturbations. This procedure has cured the certain problems of the 
unflavoured backgrounds and generated previously known smooth solutions. 
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